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various aspects of these models, namely, dynamical breaking and interrelation between ordinary 
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I. INTRODUCTION 

(N ■ 

^ I Ever since the formulation of supersymmetric quantum mechanics [1], a much deeper understanding of structure of 
. various solvable potentials in non-relativistic quantum mechanics has been gained through the ideas of supersymmetry 

■ (SUSY) (for a good review, see, e.g., [2, 3]). For instance, seemingly distinct potentials, e.g., an infinite square- well 
and a cosec^ potential, are in fact supersymmetric partners, and their respective energy levels are connected by 
supercharges. The methods of SUSY have also been applied to obtaining exact spectra of quantum systems with multi- 
component wave functions, such as Pauli and Dirac equations. It is based on the fact that, for some electromagnetic 
field configurations, Pauli and Dirac equations possess intrinsic SUSY structure, and different components of a wave 

^ function are related by some supercharges [4, 5] . More recently, the method of SUSY has been employed to construct 
quasi-exactly solvable potentials in the Pauli and Dirac equations [6, 7]. For quasi-exactly solvable systems, it is only 
O i' possible to determine algebraically a part of the spectrum but not the whole spectrum [8-12]. 

D , Recently, the concepts of SUSY and quasi-exact solvability were combined within the framework of so-called A/'-fold 
' supersymmetry [13], which is a natural generalization [14] of ordinary supersymmetric quantum mechanics. Among 
^ , various nonlinear extensions of ordinary SUSY such as parasupersymmetry [15-17], fractional supersymmetry [18], 
and so on, A/'-fold SUSY is characterized by the fact that anti-commutators of fermionic operators are polynomials 
. , of degree (at most) Af in bosonic operators. It has been proved in a generic way that A/'-fold SUSY is essentially 

■ equivalent to (weak) quasi-solvability (this latter term, in contrast to quasi-exact solvability, is used to include the 
case where the system admits non-normalizable solutions in closed form) [13]. Up to now, three different families 
of AT-fold supersymmetric systems have been found for arbitrary finite integer Af, namely, type A [19, 20], type B 
[21], and type C [22], which have correspondence with the classification of second-order linear differential operators 
preserving a monomial-type vector space [23]. Of the three types, type A and type C under a specific condition have 
been completely classified [20, 22]. 

While there are still a lot to be done in the mathematical developments of A/'-fold SUSY, it is also interesting that 
one looks for physical systems which possess such generalized SUSY. In view of the fact that so far all the A/'-fold 
SUSY potentials are only one-dimensional, it is natural that one should look for physical models which are effectively 
one-dimensional. Experience gained in the work in [6, 7] suggests that Pauli and Dirac equations are good candidates 
to start with. In this respect we note that the authors of [24] found that if the Pauli equation is generalized such that 
the gyromagnetic ratio g — 2 of the electron is changed to some unphysical values g = 2J\f {Af > 2), then for certain 
forms of magnetic fields, the generalized Pauli equation could possess type A A/'-fold SUSY. The result is interesting 
in a purely mathematical view point, but unfortunately it would not describe any existing physical systems. 
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Therefore, we would like to extend the realistie systems in [6, 7] to include A/'-fold SUSY. Since the Pauli and 
Dirac equations considered there all possess ordinary SUSY, it is therefore natural that we look for an ordinary su- 
persymmetric system which has jV-fold supersymmetry as well, as a starting point of the aforementioned purpose. 
This naturally led us further to consider a more general situation in which a system has simultaneous A/'-fold super- 
symmetry with two different values of jV. In our previous paper [25], we have succeeded in constructing a family of 
such systems, which we have called type A (A/i, A/'2)-fold SUSY. These systems possess simultaneous type A Mi- and 
7V2-fold SUSYs. Hence, the aim mentioned at the beginning of this paragraph is equivalent to finding Pauli and Dirac 
systems with type A (A/", l)-fold SUSY. 

In this paper, we shall present a detailed study of Schrodinger, Pauli, and Dirac systems which possess ordinary and 
type A SUSYs simultaneously. It is amazing to realize that some well known shape-invariant potentials, namely, the 
infinite square- well and the cosec^ potentials with specific strengths of the coupling constant, are indeed potentials 
of such kind. We also find an elliptic generalization of these potentials which breaks shape invariance but preserves 
the structure of type A (A/", l)-fold SUSY. We study various aspects of these models such as dynamical breaking and 
interrelation between ordinary and A/'-fold supersymmetries, shape invariance, quasi-solvability, and an associated 
algebra which is composed of one bosonic and four fermionic operators and dubbed type A (A/", l)-fold superalgebra. 

The organization of the paper is as follows. In Sect. II some of the main ideas of A/'-fold SUSY are reviewed. We 
then discuss in Sect. Ill the idea of simultaneous type A (A/', l)-fold SUSY. Section IV gives a detailed discussion of 
trigonometric Rosen-Morse type potentials in Schrodinger equation. In particular, we show that shape invariance in 
the case enables us to establish a relation between ordinary and A/^-fold supercharges, and to derive type A (A/', l)-fold 
superalgebra in closed form for arbitrary A/". Similar discussion of elliptic potentials is given in Sect. V. Extensions 
to the Pauli and the Dirac equations are presented in Sect. VI. Section VII concludes the paper. In the Appendix A, 
we list the complete classification of real type A (AA, l)-fold SUSY. Appendix B gives a summary of the generalized 
Bender Dunne polynomials (GBDPs), and present some explicit forms of them for the shape-invariant and elliptic 
potentials investigated in Sections IV and V. 



II. A/'-FOLD SUPERSYMMETRY 
A. Definition of A/^-foId Supersymmetry 

We begin with defining A/'-fold supersymmetry in one- dimensional quantum mechanics. Consider a matrix Hamil- 
tonian Hjv" given by 

S-)^ (2-1) 



where are ordinary scalar Hamiltonians 



Hj^ 



H^ = -p' + V^{q), (2.2) 



with p = —id/dq. A/^-fold supercharges are defined by. 



PZr \ ^+ / 



where the components are defined by, 

= Pm, P^ = Plf, (2.4) 

in terms of a monic A/'th-order linear differential operator Pa/ of the form, 

d-^ d-^-i d 

" d^ ^"'^"^^^^ + • • • + wi(g) — + w;o(9). (2.5) 

In Eq. (2.4), the superscript t denotes the formal transposition of operators in a linear space of functions of q defined 
by {d/dqY = —d/dq. Note that when all the functions Wk{q) appearing in Eq. (2.5) are real- valued, the operator P^ 
defined by Eq. (2.4) is identical with the formal adjoint of Pj^j-: P^ = P^. 
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A system (2.1) is said to be jV-fold supersymmetric with respect to Qj^ if it commutes with them: 

[Qj^,H^] = [Q+-,HAr] =0. (2.6) 
In components, the commutation relations (2.6) are expressed as the following intertwining relations: 

P^H^ - H+Pj^ = 0, P]^H+ - H^P+ = 0. (2.7) 

Therefore, the relations in (2.7) give the conditions for the system ^^_^f to be AA-fold supersymmetric. Note that the 
Hamiltonians (2.2) are always symmetric under the formal transposition even when they are not Hermitian, and thus 
each of the relations in Eq. (2.7) actually implies the other. 

It was investigated in Ref. [13] that the AT-fold supersymmetric models defined above have several significant 
properties similar to those of the ordinary supersymmetric models. In the following, we shall summarize some of the 
most important aspects of A/"- fold supersymmetry. 



B. Quasi-solvability, mother Hamiltonian, and generalized Witten index 

One of the most important aspects of A/'-fold supersymmetry is that the component Hamiltonians H~ and are 
always weakly quasi-solvable with respect to the operators PJ^ and P^, respectively [13, 20]. That is, leave the 
kernels of P^ invariant: 

H^V^CV^, V^ = kerP±. (2.8) 

As a consequence, we can in principle diagonalize algebraically the Hamiltonians H"^ in the finite Ay-dimensional 
vector spaces V^, which are thus called solvable sectors of H"^. If the space Vj^^~^ is a subspace of a Hilbert space 
L'^ on which the Hamiltonian is defined, the elements of V^'"~^ provide a part of the exact eigenfunctions and 

thus is called quasi- exactly solvable. 

In ordinary SUSY, the anti-commutator of the supercharges corresponds to the original Hamiltonian H_a/. However, 
it is not generally the case with AA-fold SUSY. This is because {Qj^, Q^v^} is now a iAAth-order differential operator. 
The half of the anti-commutator of Qj^ and Qjy)- is called mother Hamiltonian and is denoted by H^: 

An immediate consequence of the above definition is that the mother Hamiltonian always commutes with the AA-fold 
supercharges, that is, it is AA-fold supersymmetric: 

[Qj^,WA^] = [Qj^,HA^] =0. (2.10) 

Furthermore, if the original Hamiltonian Htv' is AT-fold supersymmetric, the mother Hamiltonian also commutes with 
Hjv" due to the relation (2.6): 

[Ha^,Ha^]=0. (2.11) 

We now come to discuss the characteristics of the spectrum of an AT-fold supersymmetric system. As with ordinary 
SUSY, the bosonic (lower) and the fermionic (upper) states of such a system are in one-to-one correspondence provided 
that these states arc cigcnstatcs of the mother Hamiltonian with non-zero eigenvalues. This can be seen as follows. 
Consider a normalized bosonic state $~ satisfying 

=£;-$-. (2.12) 
Since Hj^r commutes with H^v", $~ is also an eigenstate of Hj^: 

il^ = ^ il^ ■ (2-13) 



If £ is non-zero, then there exists the following normalized state 



(2.14) 



4 



Prom Eq. (2.7), we have 

Hj;^^+ = E-^+, (2.15) 

which shows that is also an eigenstate of Hj^ with the same eigenvalue E~ . Furthermore, this state is also the 
eigenstate of the mother Hamiltonian with the same £: 

(X) = ^ (^o") ' (2-16) 



since H^f commutes with Qj^ and 



(2.17) 



$- = ^=P^^$+. (2.18) 



; ^/2f V*" 

Similarly, one can show that a bosonic state is obtainable from a fermionic one at each energy level unless £ = 0: 

1 

^£ 

One can generalize the Witten index of ordinary SUSY to jV-fold SUSY by 

Tr(-l)^ = dimker Qj^ - dimkerQj^, (2.19) 

where ker are subspaces of a Hilbert space on which the super-Hamiltonian Hjv" is naturally defined. We note 

that only the physical states with energy £ = oi the mother Hamiltonian Tix contribute to the index. The index 
thus takes an integer value as the number of these zero energy states is finite {2Af at most). When the generalized 
Witten index is non-zero, AA-fold supersymmetry is not broken dynamically. 



C. Type A A/^-fold supersymmetry 

As previously mentioned, at present three different families of AA-fold supersymmetric systems have been found for 

arbitrary finite integer J\f. Wo shall be concerned only with A/'-fold supersymmetry of type A. 
The component of the type A A/'-fold supercharge is defined by 

P.. = + Wi,) - ^E(,)) (A + wi,) - ^Ei,)) X . . . 

• • • X ( A + w{q) + ^£^(g)) + Wiq) + ^^m) ■ (2.20) 

According to Ref . [20] , the necessary and sufficient condition for type A A/'-fold supersymmetry is the following: 

V^il) = Iwiqf - ^f^{2E'{q) - Eiqf) ± ^W'iq) - R, (2.21) 

where is a constant, and the functions W{q) and E{q) satisfy 

d , A d / d 

E a\ — — 



A_2i.(,))(A_^(,))iL(^ + ^(,))^(,) = for Af>S. (2.23) 

Construction of potentials possessing type A A/'-fold supersymmetry has been investigated by analytic calculations of 
the intertwining relations, and by s[(2) construction based on quasi-solvability. We refer the readers to [20] for the 
complete classification of the potentials. 

The solvable sectors of the type A A/'-fold supersymmetric Hamiltonians, which are the vector spaces preserved 
by Hj^, are given by the kernel of the type A A/^-fold supercharges: 



AT 



= kerP± = e-^^(l, z,..., z^-'), (2.24) 
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where the gauge functions VV^ are given by 

W^ = ^^^ I dqET I dqW. (2.25) 

Finally, wc note here that, when A/" = 1, 1-fold SUSY is just ordinary SUSY, and PJ^ and are simply ordinary 
supercharges A~ = + W and ^+ = + respectively. 

III. SIMULTANEOUS ORDINARY AND Af-FOLD SUSYS 

Since various quantum systems of physical interest possess ordinary SUSY, an interesting question is therefore 
whether some of these also possess AA-fold SUSY. This question has led us to investigate a more general problem, 
namely, to classify systems which possess type A A/'- fold SUSY with two different values of Af [25] . We have called 
such systems as type A (A/i, A/'2)-fold supersymmetric. 

A system H is said to have {M\,M2)-fold supersymmetry if it commutes with two different A/i-fold supercharges 
(i = 1,2) simultaneously, namely, 

[Q^r^H] = [Q(^f,H]=0. (3.1) 

Without loss of generality, we can assume A/i > Af2- In this case, it is evident that the components of the system 
H preserve two vector spaces Vy^'^ = kev pj^^^ and V^^^ = ker pj^^^ separately, where Pj^''.^ are components of 

Qjl^f" [i = 1,2). Hence, the solvable sectors V^^^j^^ of (A/i, A/'2)-fold supersymmetric Hamiltonians H"^ are generally 
given by 

VlA.,=V^)^UV;^f. (3.2) 

Classification of type A (A/i, A2)-fold supersymmetry has been presented in [25]. Here we shall focus on type A 
(A/', l)-fold supersymmetric system, i.e., system possessing simultaneous ordinary and A/'-fold SUSYs. For such a 
system, its potential must satisfy 

V^{q) = \w{qf - ^^^{2E'{q) - E{qf) ± {q) - R 

= \{Wu{qf±W'j^{q)), (3.3) 

where i? is a constant, and W^f is the first derivative of a superpotential of ordinary SUSY, with the subscript J\f 
signifying that it is also related to TV-fold SUSY. 
From Eq. (2.21) we immediately have 

MW = W'j^, (3.4) 
- ^ {2E' -E^)-2R = W^. (3.5) 

The first condition (3.3) is easily integrated as 

Wj^ = AfW + C, (3.6) 
where C is a constant. Substituting Eq. (3.6) into Eq. (3.5), we obtain 

(7^2 -1)W^ + 2MCW + ^^^~^ {2E' - E^) + C^ + 2R = 0. (3.7) 

To investigate Eq. (3.7), it is more convenient to gauge-transform the type A Hamiltonians to [20] 



F± = e-w^ 



e^^, (3.8) 
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where the gauge potentials arc given by Eq. (2.25), and A{z) and Q{z) are polynomials of at most fourth- and 
second- degree, respectively, and related to E{q) and W{q) by 

^z) = ^{z'f = aAZ^ + azz^ + 022:^ + axz + ao, (3.9) 
A!{^z) = z" = Ez', (3.10) 
Q(^z) = -Wz' = b2z'^ + biz + bo. (3.11) 

With the aid of these relations, we can show that the condition (3.7) is satisfied if 
where 

H[A]= AA" -^{A'f, (3.13) 

is an algebraic covariant called the Hessian of A [26, 27]. The solvable sectors of the type A (jV, l)-fold 
supersymmetric Hamiltonians are 

V^,i=V^UV± (3.14) 
where are given by Eq. (2.24), and Vi consists of the states 

cx e=^^« exp(±N' J dqW^. (3.15) 

We note here that the condition (3.12) is expressed as an algebraically covariant form under the projective trans- 
formations GL{2,K) (Jf = M or C) on A{z) and Q{z) (cf. Refs.[ll, 20]): 

A{z)^A{z)=/^-\^z + 5fA{^^y (3.16) 

Q{z)^Q{z) = i:^-\^z + 5fQ(^^\ (3.17) 

where a, [3,^,5 £ K and A = a(5 — /37 7^ 0. We can easily show that the polynomial A{z) of at most fourth-degree 
is transformed to one of the five canonical forms listed in Table I by the GL{2, C) transformation. The complete 
exhibition of each model for more general type A (jVi, j\/'2)-fold SUSY in this complex classification scheme is found 
in Ref. [25]. 



Case Canonical Form H[A] 



I 


1/2 





II 


22 


-3 


III 






IV 


2ix{z^ - 1) 


-V(2'' + 2) 


V 


22;3 - 322/2 - 53/2 


-32^ - 2,g2Z^l2 - 6g3Z - 3firi/16 



TABLE I: Canonical forms of ^(2) and their corresponding i?[A]. The parameters H, g2, gs € C satisfy At and 32 — 27^1 ^ 0. 

However, for physical applications where the Hamiltonians should bo real and Hermitian, it is more convenient and 
suitable to express potentials in terms of real functions. For this purpose, we should use the real classification scheme 
under the GL(2,R) transformation. In Appendix A, wc list all the inequivalent models of type A {Af, l)-fold SUSY 
in terms of real functions. Among these models, we find that only the Rosen-Morse type and some of the elliptic 
models can have physical interest, which we shall discuss in the next two sections. 
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IV. SHAPE-INVARIANT POTENTIAL 

In this section, we concentrate on the trigonometric version of Case IVa, Eqs. (A8) (AlO). The hyperbohc version 
is not of physical interest, as any clement of the solvable sectors and ijj^ is non-normalizablc and thus unphysical. 
For convenience, we shall hereafter adopt a unit system in which h = e = 2m = c = 1, so that the results obtained here 
for the Schrodinger equation can be readily carried over to the Pauli and Dirac equations. As such, the Hamiltonians 
and potentials in this and the next two sections differ from those in the previous sections and in Appendix 
A by a multiplicative factor 2. 

The trigonometric version of Case IVa can be obtained by applying a scale transformation Eqs. (Al) and (A2) 
with v < to Eqs. (A8)-(A10). We take v = —1/4 without loss of generality. With this choice, the functions which 
characterize the AA-fold {Af > 2) and the ordinary supercharges are given by (we hereafter change the variable q to x) 

1 Af 
E{x) = cot a;, W{x) = — - cot a;, Wj^{x) = — — cot a;. (4.1) 

The corresponding Hamiltonians read 

= + ^^^^^ = i [-^^-^ ^ x-M']. (4.2) 

These potentials are periodic in x, but tend to positive infinity at x = nir, n = 0, ±1, ±2 . . .. Hence, the particle is 
confined only in any one of the infinite cosec^ wells with width tt in x. Without loss of generality, we shall take the 
physical domain of interest to be < x < tt. The solvable sectors are given by 

= (sinx)~^ (l,cosa;, . . . , (cosx)"^""^), ipoOi{smx)^^. (4.3) 

In what follows, we will investigate and discuss various aspects of simultaneous ordinary and AA-fold supersymmetries, 
interplay of them, and their physical consequences. 

A. Structure of ordinciry SUSY 

First of all, due to ordinary SUSY, the Hamiltonians are factorizable as 



with the first-order operators 



It follows that 



Hj^ = ^^^j^' H+ = AJ^A+ (4.4) 



Aj,^-l + W^. (4.5) 



which corresponds to Eq. (2.7) in the case of 1-fold SUSY. If we introduce the supercharges in a matrix repre- 
sentation by 

the relations (4.4) and (4.6) are expressed as a Lie superalgebra: 

[A± , Hm] = 0, { Aj^, A+ } = H^. (4.8) 

Let us denote the normalized eigenfunctions of the Hamiltonians by tp^^^ with eigenvalues E^^^ , respectively. 
Here the subscript n = 0, 1, 2, . . . denotes the number of nodes of the wave function. It is easily proved that 
and y^, being partners of ordinary SUSY, have the same energy spectrum except for the ground state of with 



Eq^^ = 0, which has no corresponding level for V^. More explicitly, we have the following supersymmetric relations: 



E^^)^=E%\-, (4.9) 
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e^^ = «i-)-^'^^^^i^i-, (4.10) 

V'i-Tl- = (i^r^^)"'''^^4-^'+- (4.11) 

Hence Aj^ converts an eigenfunction of Hj^ into an eigenfunction of with the same energy, but with one less 
number of nodes, while docs the reverse. 

The supersymmetric element tpQ of the solvable sectors in Eq. (4.3) is normalizable in the physical domain < a; < tt 
we have considered, and corresponds to the ground state wave function ipQ^~ of V^: 

= Vo" « (sinx)^, ^j^Vr^- = 0. (4.12) 
Hence, ordinary SUSY in this system is unbroken. 



B. Two chains of shape-invariant potentials 

Next, we note that the potentials satisfy 

Vj^{x)=Vj^^2ix)+R^, Rj^=M+l. (4.13) 

This implies 

^M^M = ^M+2^N+2 + Rm, (4.14) 

or, equivalently, = HJ^^2 + ^A^- Potentials satisfying Eq. (4.13) are called shape-invariant, and are (exactly) 
solvable [2, 3]. Indeed, it turns out that belong to one of the well known shape-invariant potentials, namely, the 
Rosen-Morse-I potentials [2, 3], only with parameters being related to the integer Af. In particular, the potential Vj^ 
corresponding to A/" = 2 is simply an infinite square- well in the domain < a; < tt. Below we give a brief discussion 
of the relation between shape invariance and exact solvability, specializing to the potentials considered here. 

The wave function and energy of the ground state of Vj^_|_2 are V'o^^^^ (sin x) (cf. Eq. (4.12)) and Eq^^'^^~ = 
0, respectively. The shape invariance (4.13) indicates that the ground state wave function rp^^^+ of Vj^ is given by 

ipf^''^ = ipo^~^'^'^~ with the ground state energy Eq^'^'^ = Rj\f. On the other hand, from the supersymmetric relations 
(4.9) (4.11), the first excited state V'l'^''" of Vj^ is related to the ground state V'o^^^ of by "f/"!^^" cx A^-i/'q"^''^, 
with the same energy = Ejf^^^ = Rj^. Thus, we have ip^'^~ oc Aj^ipQ^~^'^^~ . By iterating this process, we can 

obtain the wave function of the nth excited state of Vj^ as 

and energy eigenvalues are given by 

n-l 

E^f^- = Y,R^+^k=n{N + n), n = 0,1,2,.... (4.16) 
fe=o 

With these formulas, the energy and wave functions for all V^_|_2fc are also exactly calculable in closed form. We note 
here that since Eq. (4.13) requires the integers N in to be differed by two, we therefore have two different chains 
of shape-invariant cosec^ potentials in which two neighboring potentials possess type A (TV, l)-fold SUSY, one for 
even A/^'s, and the other for odd A/''s. 

To make discussion of A/'-fold SUSY structure of transparent in the next subsection, we find it convenient to 
rewrite the symbols of the potentials in the same chain of shape invariance related by supercharges Ajj-_^2k ^^'^ ^A/'-i-2fe' 
(fc = 0, 1, 2, . . . and N + 2k> 2). Let us write 

k-l 

(a;) = V;^+2fc(^) + E = 0, 1, 2, . . . . (4.17) 

1=0 

The potentials U(") and 1/(^+2) possess type A (m, l)-fold SUSY. Now let us denote by vi"^ and ^^"^ (m = 
A/" + 2k, n,k = 0,1,2,...) the wave function and energy eigenvalue, respectively, of the nth energy level of the 
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potential V^"^\ According to Eq. (4.17), the ground state of l/^™) is also the ground state of Vj~ (with energy 
= 0), hence we have Vo"'' « (sina;)^, and i;^™' = (m^ - A/'2)/4 from Eqs. (4.12) and (4.16). Then the 
supersymmetric relations (4.4), (4.9)-(4.11) can be rewritten in terms of the eigenfunctions Vrt™^ and eigenvalues 
E^^ of the potential defined by Eq. (4.17) as follows: 

A+A-^(r^ = Eir^i;(r\ (4.18) 
Ei:i{^Eir+'\ (4.19) 

A- « (4.20) 
^+-2V'i"^«V^f"7'^- (4.21) 
The relations (4.20) and (4.21) are represented diagrammatically in Fig. 1 for odd Af. 



C. Structure of type A A/"- fold SUSY 

Now wc examine how type A A/"- fold SUSY manifests itself in the two chains of shape-invariant potentials. Es- 
sentially three issues need be addressed, namely, (1) whether AA-fold SUSY is dynamically broken or not; (2) how 
the supercharges and are related; and (3) what kind of superalgebra the five operators H^v", A^, which 
characterize the system satisfy. We shall address these questions below for the two chains of shape-invariant potentials. 



1. Dynamical J\f -fold SUSY breaking 

If AA-fold SUSY is not broken, then the ground state, and perhaps some excited states, must belong to 

Vj^ = (sinx)-^+i(l, cosx, . . . , (cosx)-^-^). (4.22) 

In our present case, the ground state of V^, given by Eq. (4.12), is arranged as 

V'^-^^ oc (sinx)^ = (sina;)-'^+i(sina;)^-^ (4.23) 

Now it is not hard to check that for odd A/', the factor (sina;)-'^"-^ belongs to (l, cos a;, . . . , (cosa;)-'^"^), while for even 
Af it docs not. 

Next, we shall consider the first excited states f/'i'^''- From the relation (4.21), we find that ijj^'' can be obtained 
by applying the supercharge A'^ on ipQ^~^'^\ i.e., 

a 4^+^) = - Y (sinx)^ 

oc (sina;)"^ cosa; = (sinx)~°^"'"^(sina;)-^~^ cosa;. (4.24) 

We easily see that Vi'^'' does not belong to VJ^ for any A/". Similarly, one can check that no higher excited states 
belong to Vj^, as they pick up higher powers in sin and cos terms as Af increases. 

In the same way, one can check that all the physical states of — do not belong to for any Af. 

Hence, from the generalized Witten index, we conclude that AA-fold SUSY is preserved barely by the ground state of 
Vj^ = vjf^ for odd Af, and is completely broken for even Af. 



2. Actions of on eigenstates 

As we have seen in Eq. (4.1), we have E = —2W = cot x in our case. This significant property enables us to express 
the A/'- fold supercharges in terms of the ordinary supercharges A^. In fact, the type A A/'- fold supercharge defined 
by Eq. (2.20) in this case becomes 
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PkT = [ ^ cot a; -; — cot X • • • — H — cot x — H — cot a; . (4.25) 

^ \Ax 2 J\dx 2 J \dx 2 J\dx 2 J ^ ' 

Comparing each factor in Pv to defined by Eq. (4.5) with WV = — A/'(cot a;)/2, one has for odd M = 2M + 1 
(M=l,2,3,...) 

^2M+l = (^l)*^^2Af+1^2M-l ■ ■ ■ ^3 ^1 ^^^s" ' ' ' ^2M-3^2M-1' (4.26a) 
-P2M+I = (~l)^^^2M-1^2M-3 ■ ■ ■ ^3 ^1 ^1^3 ' ' ' ^2M-1^2M+1' (4.26b) 

and for even = 2M (M = 1, 2, 3, . . .) 

-^2M — (~-'-)*^^2M^2M-2 " " " ^4 ^2 ^J^J^J ' ' ' ^2M-4^2M-2' (4.27a) 
-^2M = (~-'-)^^2M-2^2M-4 '"^4-^2^0 ^2 ^4 ' ' ' -^21^1-2^21^1 ■ (4.27b) 

We first discuss the action of Pj^ on the eigenstates of V^-^^ for odd ^f. As shown in the last section, for odd 
Af = 2M + 1 the ground state Vo^^^^\ which is the solution of ^2M+i^o^^'''^^ — 0' belongs to the kernel V2M+1 of 
P2M+1- namely, ^W+i^o^*'^^^'' = 0. Hence, the ground state is annihilated by both supercharges ^2M+i 

and J^2M+i- excited states tpn^^'^^ (n > 0) we have 

^2M+l''Pn'^~^^^ °^ ^2M+l(^2M-1^2M-3 ' ' ' ^3 )(^1 ^l')(^3 ' ' ' ^2M-3^2M-l) i'n'^'^^^- (4.28) 

Applying Eq. (4.21) repeatedly, we have 

{At ■ ■ ■ A+M-3^+M-i) <x V'I'Im-i. (4-29) 

since each maps the nth excited state of to the (n + l)-th state of The action of the A'^ A^^ on 

V'i+M-i is calculated using Eqs. (4.14) and (4.18) as 

AiAU^^l^_, = {A+A^ + R,) = {eI^Im-i + Ri) ^^uIm-v (4.30) 

Hence A^ A'l'ipl^^j^_-^ oc V'I+m-i- Finally, we note that all the A~'s grouped in the first bracket in the r.h.s. of 
Eq. (4.28) simply reverse the actions of the ^+'s grouped in the third bracket, and from Eq. (4.20) we have 

(^2M-l Am-3 • • • ^3") V'I'Im-I « ^i^""^'^- (4-31) 

Taken together, Eq. (4.28) impUes 

^'2M+i OC ^2M+i V-i^^^^^ OC • (4.32) 

Equations (4.28) and (4.32) can be visualized, as shown in Fig. 1, as a series of arrows of A+'s (starting with ^2M-i) 
pointing to the left, turning around with A^A^, and then pointing to the right with a series of ^~'s (ending with 

^2M+l)- 

Hence, the state ^''^ is obtainable from the excited state ijjn^^^^ {n > 0) either by the ordinary supercharge 
^2M+i' '^^ AA-fold supercharge P2m+i- ^he same argument, it can be shown that the actions of >ijjvf+i 

and P2M+1 '^''6 equivalent on the state tpl^^^^\ That all states in the pair of potential V^-'^^ and ^(^^+2) are related 
by both sets of supercharges is a manifestation of the underlying unbroken simultaneous SUSY structures inherent in 
the cosec^ potential with odd M = 2M + 1. 

For even Af = 2M, type A A/'-fold SUSY is broken. Also, the presence of Aq prevents the ^~ (m = 2, 4, . . . , A/" - 2) 
to reverse the actions of the corresponding on ijjr^\ Hence, the states in and are only related by 

the supercharges of unbroken ordinary SUSY, and not by those of broken type A A/'-fold SUSY. 
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.(3) 



a; 



a; 



(3) 
n+M 



.(5) 



(5) 

n+M- 




^(2M-1) y(2M+l) y(2M+3) 



A A 



A A 

^5 ^2M-3 



(2M-1) 
n+1 




A2M-I 



a' 

^2M-1 



(2M+1) 



4M+1 



'2M+1 



{2M+3) 
n-1 





FIG. 1: Equivalence of the supercharges P2M+1 ^-^i^l ^2m+i the state ip^ 
of sohd arrows of A^'s (starting with ^Jj^, 
with a series of A~'s (ending with j4jj,,j^j^) 



(2M+1) 



PjM+i ill Eq. (4.28) is visuahzed as a series 
pointing to the left, turning around with Af , and then pointing to the right 

Hence -PjM+iV'n ^'^"^^^ « ^sm+iV'"^^'^^^' « V'i^-[^^'- 



5. Operator relations among and for odd J\f 

The discussion in the previous subsection indicates that, for odd Af , the supercharge are proportional to 
when acting on the eigenstates of respectively. In this section, we will derive the operator relations between 
them which show that the A/'-fold supercharges are indeed factorizable into the product of and a polynomial 
in HJj or for odd N . 

Let us begin with the simplest case. For A/" = 3, the supercharges are given by (4.26a): 

P3- = -A^{A^A+) = -A^{H^ + i?i) = -{H+ + i?i)^-, (4.33) 

where we have made use of the property of shape invariance of the potentials, namely, Eq. (4.14). This shows that 
P^ is factorizable as a product of A'^ and a polynomial in Hf^ . The form of P^ is obtained by taking transposition. 
Now by repeated use of the shape invariance relation (4.14), one can check that the 5-fold supercharge P^ is given by 

F5- = Al{A^A^A+A+) = A^{H^ + R^){H^ + i?3 + i?i) 

= {H+ + R3){H+ + R3 + Ri)A^, (4.34) 

which is again factorizable as a product of A'^ and a polynomial in H^. 
We shall now prove the following formula by induction: 

Af-l / M-1 



n=0 ^ k=n ' 

Af-l , M-\ ^ 

= (-1)'" n ^2M+l + E I • ^2M+1- (4-35) 



n=0 ^ fe=Ti ' 

Suppose the latter formula holds for a given integer M. From Eq. (4.26a) and the assumption, we have 



^2A/+3 ^ ^2A/+3^2A/+1^2A'f+l 



Af-l / Af-l 



- (-l)^^"^^^2Af+3^2Af+l • n ( -^2Af+l + E ^2fe+l ) ' ^^Af+l- (4-36) 



n=0 ^ fc=n 
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Using Eqs. (4.4), (4.6), and (4.14), we can derive 

M-1 , M-1 



-^2M+3 ~ (~-'-)^'''^^2M+3^2M+1^2M+l ' ( "^SM+l + -R2fe+1 j 

M-1 . M 

= (-l)^"^^-4^M+3(-^2~M+3 + -R2M+1) • ( -f^2~M+3 + -R2fe+1 

n=0 ^ fe=n ^ 

M , M s 

= (-l)^+M2M+3 • n ( ^2M+3 + E ^2/^+1 ) ' (^.37) 



The last expression is notliing but tlic formula (4.35) with M replaced by Af + 1. Since we have already shown 
that Eq. (4.35) holds for A/" = 3, Eq. (4.33), we complete the proof of the formula (4.35) for arbitrary odd integer 
M = 2M+ 1(> 3). Equation (4.35) naturally explains the equivalent actions of -P2~m+i ^'^^ ^2m+i eigenstates, 
i.e., Eq. (4.32). 



J^. Type A (J\f, l)-fold Superalgebra 

It was proved [13, 28] that the anti-commutator of TV-fold supercharges is an jVth-degree polynomial in the super- 
Hamiltonian Hj^, the polynomial being (proportional to) the characteristic polynomial of the component Hamiltonians 
restricted in the invariant subspaces V^: 

{QJ^, Q+,} a det(F_^|^± - H^). (4.38) 

However, the direct calculation of the characteristic polynomial becomes a much harder task as Af increases. For 
type A and C A/'-fold supersymmetry, it was shown [20, 22] that the characteristic polynomials are given by the 
critical generalized Bender-Dunne polynomials [29] and are systematically calculated through a recursion relation. 
Hence, in the type A (jV, l)-fold supersymmetric case, we can also calculate them systematically using a recursion 
relation in the same way. We present some explicit forms of the critical generalized Bender-Dunne polynomials 
in Appendix B. Interestingly, for our present shape-invariant potential, four significant properties, namely, A/'-fold 
SUSY (2.7), ordinary SUSY (4.6), shape invariance (4.14), and the relation between ordinary and A/'-fold supercharges 
(4.26)-(4.27) enable us to calculate directly not only the anti-commutator of the A/'-fold supercharges but also all the 
anti-commutators among the ordinary and A/'-fold supercharges. This will be discussed in what follows. 
Let us first begin with the calculation of the anti-commutator of the A/^-fold supercharges: 

for odd J\f = 2M + 1. The calculation of the each component in this case is facilitated by the use of Eq. (4.35). In 
fact, from Eq. (4.35), its transposition and Eq. (4.4) we have 

M-1 / M-1 s M-1 / M-1 



-P2M+I-P2M+I - H ( -f^2M+l + E -'^2fc+l j • (^2M+1^2M+i) ' H ( -f^2M+l + E -^2;+l 

M-1 , M-1 s 2 

= H2M+I ■ n ( -^2~M+1 + E -^2/0+1 j = 7'2M+l(-ff2M+l)' (4-40) 
n=0 ^ k=n ' 



where 'P2M+1 is a monic polynomial of degree 2M -|- 1. In a similar way, we can easily show that -P2M+1-P2M+1 = 
ViM-yxiiitM-vx)- Therefore, the anti-commutator for odd M = 1M -\- 1 reads. 



M-1 , M-1 s z 

{Q2M+I' Q2M+1} = 'P2M+l(H2M+l) = H2M+I • Y\. ( H2M+I + E -^2fe+l j • (4-41) 

n=0 ^ fe=n ^ 

For instance, the anti-commutators for A/' = 3 and 5 are 

{Q3", } = ^3(H3) = H3(H3 + Rif, (4.42) 
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{Qs > Q^} = ^5(H5) = H5(H5 + ii3)'(H5 + R3 + Rif. (4.43) 

For even M = 2M, the supercharge (-Pa?) cannot be factorized as a product of Ajj {A'^) and a polynomial in 
as in Eq. (4.35). So the previous proof cannot be carried over straightforwardly. Let us begin with the simplest 
case, I.e., TV = 2. From Eq. (4.27), we have = -A:^A+. Hence, P} can be rewritten, with the use of the 
relation A^ A'^ = A^A^ and Eqs. (4.4) and (4.14), as 

^2^2 — ^0 {-^2-^2 — ^0 (^0 -^0 ~ Ro)-^t 

= {A^A+ - Ro)A^A+ = H^{H^ + Ro). (4.44) 

Similarly, we can calculate the component P4P4 as 

P4P4 = A2 Aq A2 A'^ A4 A2 Aq A2 = A2 j4q {A2 A2 — R2)A2 A2 Aq A2 
= A2{AoA+ ~Ro- R2)iAoA+ - Ro)AoA+A+ 
= {A2A+ - R2)A2A+A2A+{A2A+ + Ro) 

= H4 {H4 +R2f {H4 +R2+R0). (4.45) 
We shall now prove the following formula by induction: 

M-l y M-1 s 2 / M-1 s 

-P2M-P2M = H2M ■ n [^2M + ^ R2k] ■ ( -f^2M + ^ R2kj = 'P2m{H2m), (4-46) 

where 'P2M is a monic polynomial of degree 2M. Suppose the latter formula holds for a given integer M. Using the 
relations (2.7), (4.4), (4.14), and (4.27), we have 

-P2M+2-P2M+2 — ^2M^2M-^M+2-^2M-V2^2M-^M 

— ^2M^2M^^2M ~ ^2m)P2M^2M 

— ■^2mP2mP2m{^2M ~ ■R2m)^Jm- (4.47) 

Prom the assumption and the relations (4.4), (4.6), and (4.14), we obtain 

M-l / M-l X 2 / M-l 

. \ I U— D \ ,^ 

*2M 



IVl — 1 y IVl— L \ Z / ili — ± \ 

^2'M+2p2M+2 = ^2M^2M ' U ( + X] "^^fc j • ( H^m + X] "^^A; j {H2M " -R2m)^2A 

n=l ^ k=n ' ^ fe=0 ' 

M-l , M-l s 2 , M-l X 

= (^^2+m)'- n ^2^+ E • (H+M+T, R2k){H+^-R2M) 

n=l ^ k=n ' ^ fe=0 ' 

M , M s 2 / Ms 

= -f^2~M+2 ■ n ( -'^2M+2 + E ^2fe j • ( -ff2M+2 + E "^^fc j • (4.48) 
n=l ^ k=n ' ^ fc=0 ' 



The last expression is nothing but the formula (4.46) with M replaced by M + 1. Since we have already shown 
that Eq. (4.46) holds for Af = 2, Eq. (4.44), we complete the proof of the formula (4.46) for arbitrary even integer 
M = 2M(> 2). In a similar way, we can easily show that P2MR2M — ^2M(-f^2M)- Therefore, for even M = 2M we 
have, 



M-l . M-l s 2 / M-l X 

{Q2M' Qw} = ^2m(H2m) = H2M • n ( H2M + ^ ^2fe J ' ( H2M + ^ ^2fe J • (4.49) 

n=l ^ k=n ' ^ fc=0 ' 

For instance, the anti-commutators for M = 2 and 4 arc 

{Q2", Q2+} = ^2(H2) = H2(H2 + i?fl), (4.50) 

{Qr, Q4 } = ^4(H4) = H4(H4 + i?2)'(H4 + i?2 + Rq) ■ (4.51) 

It is interesting to note that the polynomial system \V]^{E)) for both odd and even M satisfies the following 
recursion relation with Vi{E) = E and Vo{E) = 1: 

Vu+2{E) = E{E + R^)Vm{E + Rm). (4.52) 
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Next, we shall consider the anti-commutators between the ordinary and jV-fold supercharges: 

For odd M = 2M + 1, we can easily calculate each component of the anti-commutators. For instance, using Eqs. (4.4) 
and (4.35), we immediately have 

Af-l . M-1 

^tu+l^^M+l — (~-'-)^^2M+1^2M+l ' n ( '^2M+1 + R^k+l 

n=0 ^ fe=n 

M , M-1 s 

= (-1)^ n ( iJ^-^+i + ^ i?2fc+i ) ^ (-l)*^7^M+l(i^2M+l), (4-54) 

where 7?.m+i is a monic polynomial of degree M -\-\ (we shall adopt the convention that a summation is zero if 

^2M+1^2M+l ~ 



the upper index is smaller than the lower one, i.e., X]n=a /« = if 6 < a). Similarly, we obtain ^2M+i^2^ 



(— l)^7?.M+i(-f?iv^_l_i) -'^w+i^Jm+i ~ (~l)^^M+i(-ff2M+i)- Therefore, the anti-commutators for odd M = 
2M -h 1 reads. 



M . M-1 V 

{AW'Q2M+l} = {A2M+1'QW} = (-l)^^M+l(H2M+l) = (-1)^ [J (h2M+1 + ^2^+0 " (^-^S) 



n— ^ k—n 



For even A/" = 2M, each component in the anti-commutators cannot be expressed as a polynomial in Hj^ or H'j^. 
However, with the aid of Eqs. (2.7), (4.4), (4.6), and (4.49) we have 

P^P+A^ = < • Vj^{H^)A-^ = . Vm{H]^), (4.56a) 
P^AJjA+Pn = P^H+P^ = ■ VMH^)- (4.56b) 

In a similar way, we can show PJ^Aj^Aj^P^ = Ajj-P^PJ^A'j^ = ■ Vm{H^). Therefore, combining the results, we 
obtain the following algebraic relation: 

{ A+ , Qj^} • {A^, Q+ } = {A^, Q+ } • {A+ , Q^} = • Vh{^m). (4.57) 

In particular, for even M = 2M we have 

M ^ M-1 ^ 2 / M-1 V 

{A^M,Q?M}-{AjM>Qw} = n(H2M+ $^i?2fe) • (h2M+ 5^ i?2fe). (4.58) 

n=l ^ k=n ' ^ fe=0 ' 

The existence of the last factor in thc> r.h.s. naturally explains why the anti-commutator of A^ and for even M 
cannot be expressed as a polynomial in H_/v'. From Eqs. (4.55) and (4.57), or directly from the expressions of P2M+1 
in Eq. (4.40) and 7?.m-i-i in Eq. (4.55), we see the following relation: 

Uu^xiEf = EV^M+i{E). (4.59) 

Finally, we summarize the complete type A {M, l)-fold superalgebra composed of one bosonic and four fermionic 
operators {Hjv', A^, Q^} for both odd and even N: 

Type A {J\f, l)-fold superalgebra: 

[A±,H^] = [Q±,H^] =0, (4.60a) 

{ A± , A± } = {A± , Q± } = {Q± , Q± } = 0, (4.60b) 

{ Aj^, A+ } = H^, {Qj^, Q+ } = P^(H^), (4.60c) 

{A^, Q^} • {A^, Q^} = • V^iU^). (4.60d) 

For odd Af = 2M + 1, the last relation can be decomposed as 

[A.^M+i^^'^M+i] = (-l)*'7^M+l(H2M+l). (4.60e) 
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(4.61) 



Using the explicit value of Rj^ in Eq. (4.13), the monic polynomials V// and TZm+i are calculated as 

E ■ Y[{E + M^ -n^y ■ (E + M^) for A/" = 2M, 

n=l 

E ■ Y[ {E + {M - n){M + n + 1))'^ forA/'=2M+l, 

. n=0 
M" 

TZm+i {E) = l[{E+{M-n){M + n + l)). (4.62) 

V. ELLIPTIC POTENTIAL 

Next, wc consider an example from Case V, i.e., the elliptic case. For physical systems, it is more relevant to 
consider this case in terms of the real Jacobi elliptic functions instead of the Weierstrass functions discussed in [25] . 
A case of physical interest, which corresponds to Case V'a in Appendix A, is given by the potentials 



^aa(^) = 2„ ± IT „„2„.„2^ - -T-2— - 1 ) . (5-1) 



E{x)= "^"f^ (2dn^x-l), Wm{x)^NW{x) . (5.2) 



where < fc < 1, and fc'^ = 1 — fc^. The functions in the jV-fold and ordinary supercharges are 

sn a; dn a; ^" ^' ' ••jyy-j .--..v-/ 2sna;dna;' 
The solvable sectors are given by 

= (sna;)-^^^ (dna:)-^^^(l,cnx, . . . , (cna;)-^-^), (5.3a) 



,4_ /sna;\T^ 

^0 oc (— ) . (5.3b) 



These potentials are periodic in x with period 2K, where K is the complete elliptic integral of the first kind satisfying 
snK = 1. As in the case of cosec^ potentials, the particle is confined in one of these periodic intervals bounded by 
infinite walls. Hence, we shall take the physical domain to be < a; < 2K without loss of generality. We mention 
here that, in the limit fc — » 0, the Jacobi functions snx, cnx, and dna; reduce to sin a;, cos a;, and 1, respectively. We 
easily see that, in this limit, the elliptic system given by Eqs. (5.1)~(5.3) exactly reduces to the trigonometric system 
(4.1)-(4.3) investigated in Sect. IV. The other limit. A: ^ 1, gives a hyperbolic case, which is not of physical interest 
as mentioned previously. In the following, we shall discuss what kinds of properties remain unchanged and what kinds 
of new features arise in comparison with the trigonometric case in Sect. IV. 



A. Structure of ordinary SUSY 



Let us begin with the structure of ordinary SUSY. The supersymmetric relations (4.4)-(4.11) are completely pre- 
served in the elliptic case with 

_ _^ N cnx 
dx 2snxdna; 

The supersymmetric element of the solvable sectors in Eq. (5.3b) remains normalizable in the physical domain 
< a; < 2K we have considered, and thus corresponds to the ground state wave function tpo^^~ of 

V'r-=^o-oc(£^)^, Aj^^r-=0. (5.5) 

Hence, ordinary SUSY is preserved in the elliptic case, too. 

Unlike the cosec^-potentials discussed in Sect. IV, however, the potentials (5.1) are not shape-invariant, as they do 
not satisfy Eq. (4.13) with a;-independent R_^f. Thus, all the remarkable properties of the trigonometric potentials 
discussed in Sect. IV B are lost. In particular, the Hamiltonians with the elliptic potentials (5.1) are only quasi-solvable 
but are not (exactly) solvable, which is also guaranteed by the fact 04 ^ (cf. Refs. [9, 10]). 
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B. Structure of A/'-fold SUSY 



Next, we shall examine the aspect of A/'-fold SUSY. The ground state wave function (5.5) can be arranged as 

oc (snx)-*+i(dnx)-'^(sna;)^-i. (5.6) 

Thus, as in the trigonometric case, we easily see that for odd J\f the ground state wave function i'^^'^ is an element 
of VJ^ given in Eq. (5.3a), but not for even M. For the excited states, we cannot proceed the investigation as in the 
trigonometric case due to the lack of shape invariance. However, we can check by using the method of Bethe ansatz 
equations [6, 7, 12] that no elements of can be normalizable physical eigenstates except for the ground state V'o^^~ 
for odd N. 

This means A/'-fold SUSY is preserved barely by the ground state of for odd M while it is completely broken 
for even M . Therefore, the dynamical aspect of both ordinary and A/'-fold SUSY breaking in the trigonometric model 
is unchanged in the elliptic model although some characteristic features in the former, such as shape invariance and 
exact solvability, are lost in the latter. 

Type A (A/", l)-fold superalgebra for arbitrary M cannot be obtained explicitly in closed form in the elliptic case. 
However, the anti-commutator of the A/'-fold supercharges in this case is still given by the critical GBDP of degree N in 
Hjv". In Appendix B, we show the recursion relation which generates the GBDPs and the first several critical GBDPs 
7r](^ for the elliptic model in Case V'a. For the anti-commutators between the ordinary and A/'-fold supercharges, we 
can derive 

{A± , Q^} • {A^, Q± } = ■ iT^^i^u), (5.7) 

using the ordinary and A/'-fold SUSY relations (2.7), (4.4), (4.6), and (B5). Hence, type A (A/". l)-fold superalgebra 
for the elliptic case is simply given by Eq. (4.60) with generalized to n^^^ (note that 7r]()'' — > as A: 0). 
Unfortunately, we cannot obtain a closed formula for the anti-commutators between and in terms of Hjv" 
for both odd and even N . However, direct calculations for smaller M indicate that the A/'-fold supercharges are 
factorizable for odd N in the elliptic case, too. Suppose, for instance, is factorizable as 

P3- = -^3 (F3- + Rc) = - A3- {A+A^ + i?e), (5.8) 

where Rc is a constant. Substituting Eqs. (2.20) and (4.5) with Wjv = MW in the latter equation, we find that E 
and W must satisfy 

2E' -E^ + UW"^ + Rc = 0, (5.9) 
E" - 2W" + (2W - E)E' + 12WW' - WE"^ + 2m^ + = 0. (5.10) 

The first condition (5.9) is identical with Eq. (3.7) for A/ = 3 with C = and R = i?c/3, and thus automatically 
satisfied. Eliminating the constant Rc in Eqs. (5.9) and (5.10), we obtain 

E" - 2W" - (4W + E)E' + UWW + 2WE^ - SW^ = 0. (5.11) 

In our elliptic case. Case V'a, E and W are given by Eq. (5.2). Using these formulas, we can easily check that the 
factorization condition (5.11) is indeed satisfied. Another evidence of the factorization comes from the form of the 
critical GBDPs for odd Af. From Eq. (B9), we easily see that all the critical GBDPs of odd degree (up to A/' = 9) 
have the following form: 

7rfif^,'\E)=EnM{Ef, (5.12) 

where TZm is a monic polynomial of degree M. The latter form is in fact a necessary condition for the factorization 
of a {2M + l)-fold supercharge as P2M+1 ~ i^^)^ '^2M+i'^'^i^2M+i)''> factorization is the case, we can derive 

in a similar way of the derivation of Eq. (4.41) 

{Q2M-1-I' QJm-mI = H2M-1-1 • '^m(H2m-m)^, (5.13) 

and thus Eq. (5.12) follows. Therefore, we conjecture that a decomposed relation for odd TV similar to Eq. (4.55) 
would hold in the elliptic case, too (see also the discussion in Section Vll). 
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VI. PAULI AND DIRAC EQUATIONS 

It is well known that, under certain field configurations, the Pauli and the Dirac equation possess ordinary SUSY 

structures [2, 5 7]. And this fact has been well made used of in the studies of exact and quasi-exact solvability of these 
equations. The results in the last two sections now allow us to construct Pauli and Dirac systems with simultaneous 
ordinary and jV^-fold SUSYs. 



A. Pauli equation 

The two-dimensional Pauli equation describes the nonrelativistic motion of a charged spin ^ particle in an external 
magnetic field in a plane. The Hamiltonian is given by 

H = {p^+ A,f + {py + Ayf + |(V X A),a,, (6.1) 

where Px and are the momentum operators, 5 = 2 is the gyromagnetic ratio, A is the vector potential of the 
electromagnetic field, and Gz is the Pauli matrix. For uniform magnetic field, = By = and B^ = B, the system is 
exactly solvable, giving the well-known Landau levels. On the other hand, the Aharonov-Casher theorem guarantees 
that for any general magnetic field Bz = B{x, y) perpendicular to the xy plane, the ground state is exactly calculable, 
owing to the existence of supersymmetry in Eq. (6.1) [2, 4]. 

Consider a magnetic field in the asymmetric gauge given by the vector potential 

Ax{x, y) = 0, Ay{x, y) = -W{x), (6.2) 

where W{x) is an arbitrary function of x. The magnetic field B has components Bx = By = Q and Bz = —W'{x). 
The Pauli Hamiltonian is then given by 

H = pl+{py- W{x))^ - W'{x)az. (6.3) 

Due to the translational invariance along y direction of the latter system, \py,H] = 0, eigenfunctions ip of H can be 
factorized as 

tp{x,y) =e'''yytp{x). (6.4) 

Here ky (—00 < ky < 00) is an eigenvalue of py, and "ipix) is a two-component function of x. The upper and lower 
components of ■ip are then governed by the Hamiltonians H- and respectively, where 

^ = ( il+ ) ' = + ^^^^^ " ^ ^^-^^ 

In this form the SUSY structure of the Pauli equation is manifest, with W{x) = W{x) — ky playing the role of the 

first derivative of a superpotential. Once the upper component of is analytically solved for a nonzero energy, the 
lower component can be obtained by applying an appropriate supercharge on the upper component and vice versa 
[2, 3]. 

The SUSY structure of the Pauli equation can be made use of in several ways. First, from the knowledge of 
shape-invariant SUSY potentials, it was found that there are four allowed forms of shape-invariant W{x) for which 
the spectrum of the Pauli equation can be algebraically written down [2]. One of the four forms gives rise to a 
uniform magnetic field. Second, from the close connection of ordinary SUSY with quasi-exact solvability, several field 
configurations were constructed giving rise to quasi-exactly solvable Pauli Hamiltonians which have underlying sl(2) 
Lie-algebraic structures [6]. 

Now, we could further construct field configurations so that the Pauli equation describes a system possessing type 
A {M, l)-fold SUSY. One simply takes W{x) to be W^{x) = -M cnx/{2snxAnx) {M > 2). The magnetic field is 
given by 



Bz{x) = -W'{x) = ( + - 1 ) . (6.6) 



Mf 1_ 

2 V sn^ X dn^ x dn^ x 

Hence, the system (6.5) describes a quantum particle moving along the y direction with the momentum ky and 
confined in one of the potential wells in the strips x G {2nK, 2(n -|- 1)K) (n G Z) formed by the magnetic field (6.6). 
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The discussions in Section V tell us that the system is qiiasi-solvable and only the ground state with zero energy 
eigenvalue can be obtained in closed form, Eq. (5.5). For the Pauli Hamiltonian, the two-component ground state 
wave function ■tpo in this case is apparently given by 



Mx) = r ° oc \d^) . (6.7) 




All the other two-component functions (j){x) = (0~(a;), (/)"'"(a;))* with (p'^ix) G V^, which are preserved by the Pauli 
Hamiltonian (6.5), are not normalizable and thus unphysical. 

In the limit k ^ 0, the system reduces to the trigonometric model in Section IV, which is shape-invariant and 
exactly solvable. For the Pauli Hamiltonian (6.5) in this limit, we immediately have 

HMx)=Ej^^-Mx), V^„(x) = rlvJ+j'^J 1 , n = 0,l,2,.... (6.8) 

where '(/'^^"'"(a;) = 0, the component wave functions ■^1"^'" and V'i'^i^ ai^e calculated using Eqs. (4.10), (4.12), and 

(4.15), and the energy eigenvalue E^!^^ is given by Eq. (4.16). It is interesting to note that if we regard the Pauli 
Hamiltonian (6.5) as a super-Hamiltonian in the matrix representation (2.1), = H, the system satisfies type A 
(JV, l)-fold superlagebra (4.60) with (note that the upper and lower components of the Hamiltonians are reversed) 

^Hv^' "H'^A "H^')- '''' 

The other cases in Appendix A can be treated in the same manner, which results in different configurations of magnetic 
field. 



B. Dirac equations 

Dirac equation which couples minimally to a stationary vector potential can be treated as in the Pauli equation 
described in the last section, since the square of the Dirac Hamiltonian in this form is proportional to the Pauli 
equation up to an additive constant. We shall therefore not repeat the discussions. Instead, we will give an example 
of Dirac system coupled non- minimally to an electric field, which is described by the Dirac-Pauli equation [7]. 

Consider the motion of a neutral fermion of spin ^ with mass m and an anomalous magnetic moment /i, in an 
external static electric field E. The fermion is described by a time-independent four-component spinor ip which obeys 
the Dirac-Pauli equation 

Hil) = £ii), H = a-p + inj3cfE + j3m, (6.10) 
where p = — iV and the Dirac matrices are given, in the standard representation, by 



«= : h p= : \ l (6-11) 





where cr are the Pauli matrices. We introduce two-component spinors x and <^ by = (X) ¥')*• Then the Dirac-Pauli 
equation (6.10) becomes 

(T ■ {p i^J.E)x — {£ + m)ip, (6.12a) 

a ■ {p + i^j,E)ip ^ {£ -m)x- (6.12b) 

When the electric field is radial, E = Er{r)r, we have a set of commuting operators {H, J^, J^, — 3/4, K} where 
J is the total angular momentum J = L + S, L is the orbital angular momentum, S — is the spin operator, and 

K is defined as K = ■ L + 1). Simultaneous eigenfunctions of these commuting operators can be written as 
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where . {6, <j)) is the spin harmonics satisfying 



1 3 









— fnjyjnij ) 


^ jrrij 









|m,|<j, (6.15) 

'« = ±(i + ^), (6.16) 

(6.17) 



With these relations, Eq. (6.12) reduces to 



^ + ^ + ^Er{r)^ f_{r) = i£ + m) /+(r), (6.18a) 
+ ^ + ^Erir)^ /+(r) = iS-m) /_(r). (6.18b) 

This shows that /_ and /+ form a one-dimensional SUSY pair, with the first derivative of a superpotential being 

W{r) = - + nEr(r). (6.19) 

r 

Now if wo choose W = Wj\f = — A/'cnr/(2snrdnr) (A^ > 2), the Dirac-Pauh equation (6.18) becomes a system 
possessing type A (A/", l)-fold SUSY in its lower and upper components, just like the situation in the Pauli equation 
discussed before. The required electric field configuration is given by 

^ , . K N cnr . , . 

liEr{r)= — . (6.20) 

r 2 sn r an r 

In this case, the system (6.18) describes a Dirac particle confined in a sphere r < 2K or between two spherical surfaces 
2nK <r <2{n + 1)K {n G N) by the electric field (6.20). From Eq. (6.18), the components /_ and /+ satisfy 

H^Mr) = A^Al.,Mr) = {S' - m')Mr), (6.21) 

where are defined as Eq. (4.5) with x replaced by r. Hence, the lowest eigenfunction ^o,k with energy eigenvalue 
£o,K = m for the given integer value of k is exactly computed as 

*,. = iK'""*^-'<''-*'), (6.22) 

where V'o^^~ is defined by Eq. (4.12). Due to jV-fold supersymmetry, the system also admits local analytic solutions 
with 

f^{r) e V_^= (snr)^^^(dnr)^^^(l,cnr,...,(cnr)-^-i). (6.23) 

As discussed in Section IV C 1, they are not normalizablc except for the configuration of ipo.^ in Eq. (6.22). 

In the limit fc — » 0, the Dirac-Pauli system (6.18) gets shape invariance discussed in Section IV B. Owing to it, 
infinite number of eigenfunctions 'ipn,K and corresponding eigenvalues £n,K {n- = 0, 1, 2, . . .) for the given integer value 
of K can be computed exactly as 



where ip'':^^^{r) = 0, ip^'^ and V'i'^i^ ^-re obtained from Eqs. (4.10), (4.12), and (4.15), and is given by 

Eq. (4.16). 

The other cases in Appendix A can be treated in the same manner, which results in difi^erent configurations of 
electric field. 
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VII. DISCUSSION AND SUMMARY 

In this paper, we have discussed the interesting issue of simultaneous ordinary and type A A/'- fold SUSYs, which 
wc have called type A {J\f, l)-fold SUSY. For realistic systems with real-valued potentials, there are essentially eight 
inequivalent type A (A/", l)-fold supersymmetric models: one is conformal (Case II), three of them are hyperbolic 
(trigonometric) including the Rosen-Morse type (Cases IVa, IVb, IV'a), and the other four cases are elliptic (Cases 
Va, Vb, V'a, V"a). Of these models, the trigonometric Rosen-Morse type and some of the elliptic models turn out 
to be of physical interest. 

We have fully investigated the trigonometric Rosen-Morse type, namely, cosec^ potentials and found that when 
they have type A (A/", l)-fold SUSY, they constitute two chains of shape-invariant potentials, one is the even J\f chain 
and the other is the odd J\f chain. The even J\f chain starts with the well-known SUSY pair of a square-well and a 
Rosen Morse potentials in [2]. For both of the chains, ordinary SUSY is always unbroken. On the other hand, A/'-fold 
SUSY is completely broken for the even Af chain while it is unbroken barely by the ground state for the odd Af chain. 
We have further shown that type A (A/", l)-fold SUSY together with shape invariance enable us to obtain the complete 
type A (A", l)-fold superalgebra, which is composed of one bosonic and four fermionic operators {H^v', A^, Q^}, in 
closed form for arbitrary Af. 

For the elliptic models, we have examined the one which reduces to the cosec^ potential in the limit fc — > 0. We have 

found that the dynamical aspects of both ordinary and A/'-fold SUSYs are completely the same as the trigonometric 
case although shape invariance and exact solvability are lost in the elliptic case. We have also shown that the structure 
of type A {M, l)-fold superalgebra is essentially preserved and that it is obtained through a recursion relation though 
we cannot obtain a closed formula for arbitrary A/". 

As physical applications, we have presented how type A (A^, l)-fold SUSY can be embedded in the Pauli and Dirac 
systems, namely, a nonrelativistic spin ^ charged particle in an external static magnetic field in 2 space dimensions and 
a Dirac particle coupled non-minimally to an external static electric field in 3 space dimensions. Some consequences 
of the symmetries have been also discussed briefly. 

Relation to some mathematical theorems in Ref. [28] also gets clearer. Applying the theorem on (in) dependence of 
supercharges to (A/", l)-fold SUSY, for instance, the necessary and sufficient condition for A/'-fold supercharges to 
be factorizable as a product of a polynomial in Hamiltonians and ordinary supercharges is 

{A^,Qa.}-{A^,Q^}. (7.1) 

In our present cases, the latter condition holds for the shape-invariant potentials in Section IV with odd Af (cf. 
Eq. (4.55)), and the A/^-fold supercharges are indeed factorized as Eq. (4.35) in the case. On the other hand, for the 
shape-invariant potentials with even Af, we have from Eqs. (4.6), (4.27), and Aq = —Aq 

and thus the condition (7.1) is not satisfied. It could be possible that the A/^-fold supercharges for even Af are 
factorizable as a product of a polynomial in the Hamiltonians and a differential operator of even order 2k (0 < k < M) 
which has no factor proportional to A^. But direct calculations for smaller even A" indicate that it would not be 
the case. In other words, A^ and for even A' would be an optimal set [28]. For the elliptic models in Section V, 
although we have not been able to (dis)prove the relation (7.1) for arbitrary odd and even Af, direct calculations for 
smaller Af also indicate that f actor iz ability of the A/'-fold supercharges would hold for arbitrary odd Af. 

To the best of our knowledge, it is the first time that both ordinary and A/'-fold SUSYs are discussed not only in 
Schrodinger equations but also in Pauli and Dirac equations. Hence, wc are expecting that further studies in this 
direction would reveal much more novel properties and rich structures in various physical systems. In what follows, 
we would like to mention a couple of interesting future issues as examples. 

Recently, the so-called PT-symmetric quantum theories [30] have attracted much attention. Applications of type 
A (Wi, A/'2)-fold SUSY to PT-symmetric theories would be straightforward since the framework of A/'-fold SUSY does 
not rely on the concepts of (pseu(lo-)Hermiticity and so on. For this aim, one should simply employ another suitable 
classification scheme different from the real classification in this paper. 

The present approaches to embed A'-fold SUSY into the Pauli and Dirac systems essentially rely on the intrinsic 
structure of ordinary SUSY. On the other hand, there exists a Dirac system which has no intrinsic SUSY structure 
but is quasi- (exactly) solvable [31]. One of the peculiar aspects of the system is that the underlying algebraic structure 
is the Lie superalgebra osp(2|2) instead of s[(2) [32]. Hence, it is interesting to clarify whether or not we are able to 
embed A/^-fold SUSY into such a system, too. 
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APPENDIX A: CLASSIFICATION OF REAL TYPE A {J\f, l)-FOLD SUPERSYMMETRIC MODELS 



In this appendix, we shall present a detailed classification of real type A (Af, l)-fold supersymmctric models. For 
this purpose, we must use the classification scheme of real polynomial of fourth-degree under the real projective 
transformation GL{2,'R). In Table 11, we show the real canonical forms of A{z) (second column) and the complex 
projective transformations which convert the real canonical form to the corresponding complex one in Table I (third 
column). 



Case 


Real Canonical Form 


Transformations 


I 


1/2 


Identity 


II 


2z 


Identity 


III 




Identity with fj, = ±.v 


IV 


±2i'{z^ - 1) 


Identity with /u = ±1^ 


IV' 


±2v{z^ + 1) 


/? = 7 = 0, (5 = ia with ji = ±u 


V 


±|(l-«2)(l_fe2^2) 


7 = Q, 5 = {Ak^ - fc'2)/3/(4 + fc'2) 


v 


±|(l-«2)(fc'2+fc2^2-) 


7 = Q, 5 = -(1 + 4fc2)/3/(l + 4fc'2) 


V" 


±^{l+Z^){l + k'^Z^) 


7 = ia, J = i(4/c'2 - A;2)/3/(4 + fe^) 



TABLE II: Real canonical forms of A{z) and complex projective transformations which convert them to the corresponding 
complex ones in Table I. The parameters v, k, k' €^ satisfy i'>0, 0<fc<l and k^ + k'^ = 1. 

We note that the choice of the canonical forms is not unique, and we select them so that the change of variable 
determined by Eq. (3.9) is as simple as possible; this is justified by the fact that the final form of the potentials does 
not depend on the choice of the canonical forms, cf. Ref. [20], Section 4. As a consequence, our canonical forms are 
different from, e.g., those in Ref [26]. 

Furthermore, we note that by Eq. (3.9) a rescaling of the parameters a^, 6^, R in the models by an overall nonzero 
constant factor u has the following effect on the change of variable: 

z{q; vai, vbi, vR) = z{\fvq; ai, hi, R). (Al) 

Prom this equation and Eqs. (2.25), (3.8), (3.10), and (3.11), we easily obtain the scaling relations: 

E{q; vat, fbi, vR) = sfvE{^q\ Oj, 6j, R), (A2a) 

W(g; va^, vh,. vR) = ^fvW{^fvq\ ai, h, R), (A2b) 

W{q; va„ vh,,vR) = W{Vl^q; a„ h, R), (A2c) 

V^{q;va„vh„vR) = vV^ {^/vq;a„h, R). (A2d) 

Therefore, we can set v = 1 hi the canonical forms III-V" without loss of generality; the models corresponding to an 
arbitrary value of v follow easily from Eqs. (Al) and (A2). 

In all the cases, the condition (3.12) can be satisfied if and only if C = 0. In Case IV, the condition (3.12) requires 
a further separate analysis depending on whether 60 = or 61 = 0. Similarly, in the real elliptic cases, namely, Cases 
V-V", we must make a further separate analysis depending on whether 62 = 60 = or 5i = 0. Some of them however 
have no real solutions of Eq. (3.12), for which we shall omit the presentation of the models. We also neglect Cases I 
and III since both cases lead to a trivial model, as has been already shown in Ref. [25] . 



1. Case II: A{z) = 2z, z{q) = 



Parameters: 



&2 = 5i = 0, 60 = 1, R = 0. 



(A3) 
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Supercharge: 



Potentials: 



Solvable sectors. 



Parameters: 



Supercharge: 



Potentials: 



Solvable sectors. 



Parameters: 



Supercharge: 



Potentials: 



Solvable sectors. 



E{q) = I, W{q) = -1. (A4) 



Vi?>±=,-^<1,,» ,^l''-">, *„±(,)o<,T*^ (A6) 

2. Case IV: A{z) = 2{z^ - 1), z{q) = cosh2g 
a. Case IVa: bo = 



7V2 - 1 

62 = 60 = 0, 61=2, R = . (A7) 



, 2cosh2o , cosh2g , 



V^'^^ = (sinh2g) ^^(l,cosh2g, . . . , (cosh2g)-'^"i), (AlOa) 

V'^(g) a (sinh2g)^T. (AlOb) 

6. Case IVh: 61 = 



AA^ - 1 

62 = 61 = 0, 60 = 2, R= — - — . (All) 



^'^ 2sinh^2g 2sinh% ^ ^ 



V]/'^ = (smh29)~^(tanhg)T^(l,cosh2g,...,(cosli2g)-^-i), (A14a) 
V'^(g) (X (tanhg)=F^. (A14b) 



3. Case IV': A{z) = 2{z^ + 1), z{q) = sinh2g 



Parameters: 



Supercharge: 



Potentials: 



Solvable sectors: 



a. Case IV'a: 60 = 



62 = 60 = 0, 61 = 2, R- 



2sinh2g sinh2g 



cosh 25 ' 



cosh 2q ' 



2cosh^2g 2 



V^^=^ = (cosh2g)-^^^(l,sinh2g, ...,(sinh2g)-^-i), 
V'^(g) oc (cosh 2g) "Ft. 

h. Case IV'b: 61 = 
In this case, there are no real solutions of Eq. (3.12) for 60. 

4. Case V: A{z) = i(l - a^)(l - fe^z^), z{q) = snq 
a. Case Va: 62 = &o = 



Parameters: 



Supercharge: 



Potentials: 



Solvable sectors: 



62=60=0, &i=^, ^=^^^42"^^+^'^- 



cn g dn g 



2 cngdng 



8 cn^ q dv? q 4 \ dn^ q cn^ g dn^ q 



+ 



"Kv =(cng) = (dng) 2 (l,sng,...,(sng)-' 
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b. Case Vb: 61 = 



In this case, we have two sets of real solutions of Eq. (3.12), namely, 62 = |(1 ± k) and bo = —-^{k ± 1). Here we 
only show the result corresponding to the upper-sign solutions. 

Parameters: 



62 



Supercharge: 



Potentials: 



Solvable sectors: 



k{l + k) 



bi = 0, 60 



sng(2dn^g-fc^^) 
cnqdnq 



1 + k 



W{q) 



AT'^ - 1 

R = -^^——{i + Qk + e). 
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fc(l + k)sn^q-k-l 



2 cn g dn g 



•^'^ ~ 8cn2gdn% ^ 8 V2cn2gdn% 2dn^g 



V^^^ = (cn q dn q) (1 + fc sn g)=Fi (1 + sn g)=Fi (l, sn 9, . . . , (sn 9)^-1 ), 



% (9) <^ 



cn qr dn 



(1 + fcsn5)(l + snq) 



(A23) 
(A24) 

(A25) 

(A26a) 
(A26b) 



5. Case V: A(z) = \(l - z''){k'^ - k^z^), z(q) = cnq 



a. Case V'a: 62 = &o = 



Parameters: 



Supercharge: 



Potentials: 



Solvable sectors: 



62=60 = 0, b^ = ~, R=^^[l-l{k'^-k% 



^^^^^ cng(2dn%-l) ^ = 



^^1 



sn dn 5 



jV^cn^q Af f 1 



2 sn g dn g 



i sn^ g dn g 4 \ sn-^ g dn g dn g 



- 1 



Vif)^ = (sng)-^(dng)-^(l,cng,...,(cng)^-i) 



sn g 
dng 



(A27) 
(A28) 

(A29) 

(ASOa) 
(A30b) 



6. Case Vb: 61 = 



In this case, there are no real solutions of Eq. (3.12) for 62 and 60. 
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6. Case V": A{z) = §(1 + z^){l + k'^z^), z{q) =tnq 



a. Case V"a: 62 = 6o = 



Parameters: 



Supercharge: 



Potentials: 



Solvable sectors: 



b2 = bo = 0, &i = ^, R = -^^^^^^{l + k'^). 



2 



^, , sno(A:^ cn^ g + 2A:'^) , /e^sngcng 

E{q) = ^ W{q) = -- 



cn g dn g 



2dng 



dn^g 



dn^g 



dn^g 



V^^'± = (dng)-^(cngy^-i(l, tng. . . . , (tng)^"!), 
V'^(g) (X (dng)±^. 



(A31) 



(A32) 



(A33) 



(A34a) 
(A34b) 



b. Case V"b: 61 = 

In this case, there are no real solutions of Eq. (3.12) for 62 and 60. 

APPENDIX B: GENERALIZED BENDER-DUNNE POLYNOMIALS 

In this appendix, we summarize the generalized Bender-Dunne polynomials (GBDPs) and present some explicit 
forms of them for the shape-invariant and elliptic potentials investigated in Sections IV and V. They are a natural 
generalization of the ones first constructed for the quasi-exactly solvable sextic anharmonic oscillator by Bender and 
Dunne [29] to all the type A A/'-fold supcrsymmetric models. For more details, sec Rcf. [20], Section 6.2. In the unit 

system adopted in Sections IV and V, the GBDPs 7rjf^'(£') are characterized by the following four-term recursion 
relation: 



(i?) = {E + 2A[fl)4^1(i^) - 4n(n - M)B^^K\:'\{E) 
+ 8n{n - l)(n - M){n -M- 1)C^^t:'^^}^{E), 



where , , and Clt'^ are given by 



n{n-J\f+l) + 
« - y jaiTbo 



{Af-l){J\f-2) 



6 



02 T 1 " 7: — ]bi+ R, 



^ 1 



n — I ai =F Oo 



(Bl) 

(B2) 
(B3) 

(B4) 



and the parameters Oj and bi correspond to the coefficients of the polynomials A{z) in Eq. (3.9) and Q{z) in Eq. (3.11), 
respectively. The degree J\f polynomial (E) in each polynomial system {Trl;^'^' {E)}^^„ is called Afth. critical GBDP. 
It was shown in Ref. [20] that the anti-commutator of the type A A/'-fold supercharges is (proportional to) the AAth 
critical GBDP in the type A Hamiltonian: 



(B5) 
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To obtain the recursion relation (Bl) for a specific case, we must first make a GL{2, C) transformation so that 
the coefRcient Uq in A{z) vanishes. For the elhptic potential in Section V, Case V'a, the canonical form is A{z) = 
5(1 - z'^){k'^ - k'^z'^). Thus, we transform A{z) as 

A{z) = A-^{z + lfA(^-^ = \z[z^ + 2{k'^ -k'')z + l]. (B6) 

In our Case V'a, Q{z) = —z/2 from Eq. (A27) and it is transformed simultaneously as 

Q{z) = A-\z + J^) = - !)• (B7) 

Hence, the values of the parameters which should be substituted into the recursion relation (Bl) are as follows: 

as = ^, (12 = fc'^ - k"^, ai = ^, 04 = ao = 0, 

b2 = -j, bi = 0, bo = ^, 

and R = {k'"^ — k'^){N'^ — 1)/12 (cf. Eq. (A27)). With these parameter values, the recursion relation (Bl) reduces to 

- fc2 



E+ {2n-M+l) 



4^(^) -n(n-^)(n- ^) (n- (E). (B8) 



Using this three-term recursion relation, we can calculate, in principle, the critical polynomial tt^ {E) for an arbitrary 
value of N ■ For instance, we have 

TT™ (£)=£;, (B9a) 
^\^\E)={E-e){E + k'^). (B9b) 
^f{E)=E[E + 2{k'''-k^)\', (B9c) 
Tri"*! [E) = [E^ + 2(2fc'2 - 3k^)E + 9k^] [E^ + 2{3k'^ - 2k^)E + 9k'^] , (B9d) 
Trf\E) =E[E^ + 10(fc'2 - e)E + 8{3k'^ - 2k'^k^ + 3k^)]\ (B9e) 
ir^iE) = [E^ + (16fc'2 - 19k^)E^ + {6Ak'^ - SOk'^k^ + n5k^)E - 225k^] 

X [E^ + (19fc'2 - 16k^)E^ + ill5k'^ - 80fc'2fc2 + 6Ak^)E + 225k'^] , (B9f) 
tt^^\e) =E[E^ + 28(fc'2 - k^)E^ + 28{9k"^ - lOfc'^A:^ + 9k^)E 

+ UA{k'^ - e){5k'^ + 2fc'2fc2 + 5k^)]\ (B9g) 
TT^iE) = [E^ + 4(10fc'2 - 11/^2)^3 _^ 2(264fc'4 - 376k'^k^ + 3A7k^)E^ 

+ 4(576fc'^ - 672k'^k^ + 82Qk'^k'^ - nbbk'^)E + 11025fc^] 

X [E^ + 4(llfc'2 - lGk'^)E^ + 2(347A:"' - 'ilQk'^k^ + 264A:'')£;2 

+ 4(1155fc'*^ - 826fc'*fc^ + 672fc'^fc^ - f>7&k^)E + 11025fc'^] , (B9h) 
T:f{E) =E[E^ + 60(fc'2 - A;2)i;3 + l2(109fc'^ - lA&k'^k'^ + imk^)E'^ 
+ 16(A;'^ - fc^)(761fc'* - llSfc'^A;^ + 761fc*)i; 

+ 1152(35A;'« - 2Qk'^e + ISk'^k^ - 20k'^k^ + 35fc^)] ^ (B9i) 



In the limit fc — > (fc' ^ 1), the trigonometric case discussed in Sect. IV is recovered, and each of the above critical 



GBDP ■k{P{E) exactly reduces to the polynomial Vj^{E) in Eq. (4.61). It is interesting to note that all the above 



ir^^ have the following forms: 



4'i^::\E)=En'^\E;k^k'r, (BIO) 
nf;^\E)=TZ^^\E;k^k'')TZ^^\E;-k'^-k% (Bll) 
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where TV-^i and TZ ^^ are monic polynomials of degree M in the first variable, though we have not been able to prove 
them for arbitrary M . 
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